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Markierungshügel und 
Zeitreihenanalyse - oder: Wann beginnt 
ein Biberjahr (Castor fiber L.)? 


von B. Klenner-Fringes 


Abstract 

Scent marking and time series 
analysis or: When is the beginning of 
a beaver-year (Castor fiber L.)? 

Time series analysis is a Statistical method 
that can help to reveal rhythmic seasonal 
changes in behaviour. Düring a research pro- 
ject the question arose whether there is a 
year-to-year change in the utilization of 
riverbank structures by Castor fiber albicus 
MATSCHIE 1907. In Northwest Germany 
beavers scent-mark full-year with varying in- 
tensity. So scent marking as non-invasive 
method was used to estimate the beginning 
of a “beaver-year”. As the original data did 
not show a clear seasonal periodicity, expo- 
nential smoothing was used to visualise this 
rhythmical changes. 

The smoothed series reveales a clear 12- 
month-period, and in January, the number 
of scent marks exceeds the five-year-median 
and even reaches the upper 75 %-level for 
the first time. 

Knowing that mating Starts in January 
and together with the fact that since the 
end of December of the first year after the 
releasing the territories did not further en- 
large, the results of the time series analysis 
are an additional indication that January 
can be regarded as the beginning of the an- 
nual cycle in Northwest Germany. 

Key words: beaver, Castor fiber , scent 
marking, time series analysis 

Zusammenfassung 

Die Zeitreihenanalyse ist eine statisti¬ 
sche Methode, mit deren Hilfe rhythmische, 


saisonale Verhaltensänderungen verdeut¬ 
licht werden können. Im Rahmen eines Frei¬ 
landforschungsprojektes stellte sich die Fra¬ 
ge, ob sich die Nutzung von Uferstrukturen 
durch Castor fiber albicus in den einzelnen 
Jahren unterscheidet. Da Biber in Nordwest¬ 
deutschland mit unterschiedlicher Intensität 
das ganze Jahr hindurch markieren, wurde 
die Quantifizierung des Markierverhaltens 
als nicht-invasive Methode benutzt, um den 
Beginn eines „Biber-Jahres“ festzulegen. 

Da die über einen Zeitraum von fünf Jah¬ 
ren nach der Wiederansiedlung gesammelten 
Originaldaten keine deutliche Periodizität 
erkennen ließen, wurden die Daten exponen¬ 
tiell geglättet, um die rhythmischen Schwan¬ 
kungen besser sichtbar werden zu lassen. Es 
ergab sich eine 12-monatige Periodizität, wo¬ 
bei im Januar die Zahl der Markierungshügel 
den Median des gesamten Untersuchungszei¬ 
traumes überschritten und erstmals die Ober¬ 
grenze des oberen Quartils erreichte. Zusam¬ 
men mit der Tatsache, dass im Januar die 
Paarungszeit beginnt und der Beobachtung, 
dass Ende Dezember im Jahr der Ansiedlung 
die Reviere ihre maximale Ausdehnung er¬ 
reichten, können die Ergebnisse der Zeitrei¬ 
henanalyse als zusätzliche Hinweise dafür an¬ 
gesehen werden, dass der Januar als Startmo¬ 
nat für das „Biber-Jahr“ in Nordwestdeutsch' 
land angesehen werden kann. 

Einleitung 

Im Rahmen eines Freilandforschungs¬ 
projektes der Arbeitsgruppe Ethologie der 
Universität Osnabrück wurde vor einigen 
Jahren an einem Fluss im nordwestlichen 
Niedersachsen eine kleine Elbebiber-Popu¬ 
lation etabliert. 
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Ein Ziel Jer Untersuchungen war die 
Beantwortung der Frage, welche Bedeutung 
die Ufemmkturen und Ressourcen eines 
vom Menschen *tark beeinflussten größeren 
FÜeßgew&ssm in einer intensiv genutzten 
Kulturlandschalt für ein se ml aqua tische s 
Saugetier wie z, R. den Biber haben, in wel¬ 
cher Weise $v)wohl natumahe als auch an¬ 
thropogene Strukturen und Ressourcen ge¬ 
nutzt werden (Ahk 1 und 2). 

Die Auswertung der über einen Zeitraum 
von tunt Jahren gesammelten Daten zeigte 
hoehsignifikante Korrelationen zwischen 
spezifischen Verhaltenskomplexen und be- 
stimmten Uferstrukturen. Der lange Unter- 
suchungszeitnmm legte es nahe» zu überprü- 
ten, ob es eine Änderung der Präferenzen im 
Laufe der Zeit gab. Es war deshalb notwen¬ 
dig» Zeitintervalle festzulegen. Am sinnvoll' 
scen schien im vorliegenden Fall ein annueb 
les Intervall zu sein, da innerhalb dieses Zei¬ 
traumes nahezu alle Verhaltensweisen und 
physiologischen Prozesse des Bibers ublau- 
ten. Um die Daten nicht einfach nach Ku- 
lenderfahren zusammenzulassen, sondern 
nach einem von Jen Tieren bestimmten Jah- 
resrhvthmus. wurde nach einer Methode ge¬ 
sucht, ein solches „Biber-Jahr“ festzulegen. 

Im folgenden soll nicht auf die Ergeb¬ 
nisse der gesamten Untersuchung eingegan- 
gen werden. Es soll lediglich eine Möglich¬ 
keit aufgezeigt werden, anhand von Biofak¬ 
ten eine annuelle Rhythmik festzulegen. 

Da cs nicht möglich w ar, bestimmte In¬ 
dividuen regelmäßig zu fangen und ihnen 
Blut abzunehmen um den Hormontiter - 
bezogen auf die entsprechenden Ge¬ 
schlecht Hormone - und damit z. B. ein über 
den Fortpftanzungsryklus definiertes Jahr 
fest zu legen, musste nach einer anderen» 
non-invasiven Möglichkeit gesucht werden, 
den Beginn eines Biberjahres festzulegen, 
etwa anhand der Zahl bestimmter Biofakte. 
Am exaktesten, Jas heißt auf ein bis zwei 


Abb, 1: Das Untersuchungsgebiet im west¬ 
lichen Niedersachsen, 

Abb. 2: Biber vor dem Bau, 

Abb. 3: Markierungshügel. 

(Alle Fotos stammen von der Verfasserin.) 
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Monat 


Abb. 4: Anzahl der 
Markierungshügel pro 
Monat. (J: Januar) 


Tage genau, lassen sich die scent marks oder 
Markierungshügel (Abb. 3), die an den 
unterschiedlichen Stellen im Biberterrito 
rium zu finden sind, zeitlich zuordnen. Es 
handelt sich hierbei um olfaktorische Mar- 
kierungen, deren Intensität im Jahresverlauf 
Schwankungen unterliegt, wie bereits 
Townsend (1953) und Hay (1958) für den 
Nordamerikanischen Biber Castor canaden - 
sis berichten. Weitere Arbeiten bestätigen 
ihre Ergebnisse, zeigen aber auch, dass es of- 
fensichtlich je nach der geografischen Lage 
der Untersuchungsgebiete Unterschiede 
gibt, und zwar sowohl für C. canadensis 
(Brenner 1967; Müller-Schwarze & 
Heckman 1980; Svendsen 1980), als auch 
für den Eurasischen Biber C. fiber (HEIDE- 
cke 1974/75; Rosell & Nolet 1997), wes- 
halb die Ergebnisse anderer Untersuchung 
gen, die in Gebieten mit abweichenden kli¬ 
matischen Verhältnissen, vor allem im Hin¬ 
blick auf die Temperatur, nicht einfach auf 
die vorliegende Untersuchung übertragen 
werden dürfen. Es war deshalb notwendig, 
für das Untersuchungsgebiet im nordwest¬ 
lichen Niedersachsen eine Analyse anhand 
der dort gewonnenen Daten durchzuführen. 

Material und Methode 


Eine Prognose sollte nicht erstellt wer¬ 
den, ebensowenig war das Ziel der Auswer¬ 
tung die Berechnung eines Trends, da für die 
Fragestellung (Wann beginnt ein Biberjahr?) 
eine mögliche Zunahme der Anzahl der 
Markierungshügel infolge der steigenden 
Zahl der Individuen und damit der Reviere 
im Laufe der Zeit nicht von Bedeutung war. 
Vielmehr kam es darauf an, saisonale, über 
den Gesamtzeitraum hin relativ konstante 
Schwankungen sichtbar zu machen und da¬ 
nach den Beginn des Biberjahres festzulegen. 

Um geringfügige, evtl, durch Witte¬ 
rungseinflüsse verursachte Schwankungen 
im saisonalen Rhythmus auszugleichen, 
wurde die Zahl der Markierungshügel für je¬ 
den Monat zusammengefasst. 

In die Auswertung gehen Daten aus 46 
Monaten des Untersuchungszeitraumes von 
fünf Jahren mit 596 Markierungshügeln ein. 
Bei den Leerwerten (s. Abb. 4) handelt es 
sich entweder um Nullwerte, d. h. in diesen 
Monaten wurden keine Markierungen ge¬ 
funden, oder aber um arbeitstechnisch sowie 
durch langanhaltende Hochwasser bedingte 
Leerwerte. 

Die statistischen Berechnungen wurden 
mit Statistica 5.0 durchgeführt. 


Als Methode für die Festlegung eines 
„Biberjahres“ bietet sich die Zeitreihenana¬ 
lyse an, mit deren Hilfe Prognosen erstellt, 
Trends aufgezeigt und - wichtig für den vor¬ 
liegenden Fall - die Formen saisonaler 
Rhythmen sichtbar gemacht werden kön¬ 
nen (StatSoft, Inc. 1996). 


Ergebnisse 

Wie Abbildung 4 zeigt , lassen die starken 
Schwankungen zunächst kaum eine Saiso¬ 
nal ität erkennen. Aus diesem Grunde wurde 
eine exponentielle Glättung der Kurve vor¬ 
genommen. Die fehlenden Monatswerte 
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Abb. S: Vergleich der beobachteten Werte 
inklusive der durch Interpolation von 
benachbarten Werten ersetzten missing 
data (dünne Linie) und der durch 
exponentielles Glätten ermittelten Kurve 
(dicke Linie). (J: Januar) 


Abb. 6: Exponentiell geglättete Kurve des 
Verlaufs der Markierhäufigkeit und statisti¬ 
sche Kennzahlen (Median, Q75 bzw. Q25: 
oberes bzw. unteres Quartil). (J: Januar) 


(Leerwerte, missing data) wurden durch 
Interpolation von benachbarten Werten er¬ 
setzt. Da rhythmische, saisonale Schwan¬ 
kungen vorausgesetzt wurden, wurde dieses 
Verfahren dem Ersetzen der missing data 
durch den Gesamt-Mittelwert vorgezogen. 
Folgende Modellparameter wurden gewählt: 

Aus den bereits erwähnten Gründen 
wurde keine Trendberechnung durchgeführt, 
da es ausschließlich um die Festlegung des 
ersten Monats eines Biberjahres ging. Als 
additive Saisonkomponente wurde ein lag 
von 12, entsprechend 12 Monaten, gewählt. 

Die Glättungsparameter a und ß wur¬ 
den durch programminteme, automatische 
Suche bestimmt. 


Als Lack-of-fit-Indikator, der die Anpas- 
sungsgüte der geglätteten Kurve (Grad der 
Fehlanpassung) angibt, wurde der am häufig¬ 
sten verwendete mittlere quadratische Feh¬ 
ler (mean squared error) voreingestellt 

Das Ergebnis des exponentiellen Glät- 
tens ist in Abbildung 5 grafisch dargestellt. 
Anders als die Originaldaten lässt die Kurve 
der geglätteten Daten die zugrundeliegende 
12monatige Rhythmik sehr gut erkennen. 
Diese 12monatige Rhythmik zeigt sich - wie 
mit unterschiedlichen lags durchgeführte 
Tests gezeigt haben - unabhängig von der 
voreingestellten Zahl der lags. 

Entsprechend den oben genannten Pa¬ 
rametern wurde für den Glättungsparameter 
a der Wert 0,005 ermittelt. Das bedeutet, 
dass der Einfluss vorhergehender Beobach¬ 
tungen langsam abnimmt; es entsteht eine 
geglättete, angepasste Kurve, die durch Stö¬ 
rungen nur wenig beeinflusst wird. Für Ö er¬ 
gibt sich der Wert 0, d. h. es gibt eine kon¬ 
stante, unveränderliche Saisonkomponente. 
Für S 0 wird der Wert 13,06, für den Lack-of- 
Fit-Indikator der Wert 0,04 errechnet. 

Der Wilcoxon-Test für gepaarte Stich¬ 
proben wurde verwendet, um die Ähnlich¬ 
keit der beobachteten und der geglätteten 
Wertekurven zu testen. Das Ergebnis (t = 
809; Z = 1,1744; p = 0,2403) zeigt, dass die 
Kurven sich nicht signifikant unterschei¬ 
den, die beiden Datengruppen also tatsäch¬ 
lich aus der gleichen Grundgesamtheit 
stammen. Dies kann als Hinweis auf die 
Anpassungsgüte der Glättungskurve angese¬ 
hen werden. 

Um den ersten Monat des Biberjahres 
festzulegen, wurde überprüft, wie sich die 
Monatswerte der geglätteten Reihe zum Ge¬ 
samtmedian dieser Datenreihe verhalten. 
Median und Perzentile wurden dabei als La¬ 
gekriterien dem Mittelwert und der Stan¬ 
dardabweichung vorgezogen, da die Werte 
der geglätteten Kurve sowohl nach dem Kol- 
mogorov-Smimov-Test (p < 0,05) als auch 
nach dem Shapiro-Wilks-Test (W = 0,9271; 
p = 0,0012) nicht normalverteilt sind. 

Die Autokorrelationsfunktion zeigt, dass 
es mit einem Korrelationskoeffizienten von 
r = 0,799 eine starke Korrelation für lag 12 
gibt. Es liegt also eine saisonale 12monatige 
Komponente vor, deren Rhythmik als gesi- 
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chert angesehen werden kann. Außerdem 
zeigt diese Funktion, dass es kein weißes 
Rauschen gibt, der Standardfehler ist klein 
und konstant; d. h. es handelt sich nicht um 
eine zufällige Zeitreihe. 

Damit ist aber noch nicht der Beginn ei' 
nes „Biberjahres“ festgelegt, sondern ledig' 
lieh eine saisonale Rhythmik aufgezeigt. 
Um hier einen Anhaltspunkt zu finden, 
wurde überprüft, wie sich die einzelnen Mo' 
natswerte der geglätteten Reihe zum Ge- 
samtmedian verhalten. 

Wie Abbildung 6 zeigt, ist der Januar- 
Wert der erste Wert, der nicht nur den Ge- 
samt'Medianwert von 12,02 übersteigt, son- 
dem der auch an die Grenze des oberen 
Quartils (Q 75 = 14,43) heranreicht. 

Berücksichtigt man, dass im Januar die 
Paarungszeit beginnt, kann dies als zusätz- 
licher Hinweis dafür angesehen werden, das 
Biberjahr mit dem Monat Januar beginnen 
zu lassen. 

Diskussion 

Der vorliegende Artikel zeigt den Ver- 
such, mit Hilfe des Markierverhaltens der 
Biber für die jahresbezogene Auswertung 
von Daten den Beginn eines „Biberjahres“ 
festzulegen. 

Die Zeitreihenanalyse ergab, bestätigt 
durch die Autokorrelationsfunktion, dass 
tatsächlich eine zwölfmonatige, also circ- 
annuale Rhythmik vorausgesetzt werden 
kann, was bei einem k-Strategen wie dem 
Biber mit seiner intensiven, etwa zwei Jahre 
dauernden Jungtierfürsorge zu erwarten war. 
Da es bei der Analyse nicht um die internen 
oder externen Zeitgeber ging, die diesen 
12 monatigen Rhythmus bestimmen, muss 
hier nicht weiter darauf eingegangen wer- 
den. Vielmehr sollte hier entschieden wer- 
den, zu welchem Zeitpunkt der Zyklus be- 
ginnen soll. 

Warum für das Untersuchungsgebiet im 
westlichen Niedersachsen der Januar als 
Startmonat ausgewählt wurde und nicht et- 
wa der Juni, in dem - nach der geglätteten 
Kurve - am intensivsten markiert zu werden 
scheint, soll im folgenden kurz begründet 
werden. 


Der Januar weist eine Zunahme der Mar- 
kierintensität auf, die das obere Quartil der 
statistischen Kennzahlen erreicht. Doch dies 
allein wird nicht als ausreichend erachtet, 
die Daten jeweils von Januar bis Dezember 
zusammenzufassen, d. h. mit dem Januar als 
erstem Monat zu beginnen. Auch die aus- 
schließliche Berücksichtigung des Biofaktes 
Markierungshügel, also nur die statistische 
Untersuchung der saisonalen Schwankun- 
gen des MarkieVerhaltens als alleinigem 
Kriterium, könnte als unzureichend angese- 
hen werden. Zwei weitere Fakten unterstüt- 
zen allerdings die Wahl des Januars: 

Gegen Ende Dezember des ersten Jahres 
der Untersuchung, 3 Monate nach dem 
Aussetzen der Biber (die Aussetzung erfolg- 
te am 12. Oktober), endete die explorative 
Phase. Die Tiere erweiterten ihren Aktions¬ 
radius danach nicht mehr (KLENNER-FRIN- 
GES 1992). 

Nach Angaben von HINZE (1961) und 
HEIDECKE (1974/75) beginnt im Januar die 
Paarungszeit. Bei Farmtieren beobachteten 
ZUROWSKI & Doboszynska (1975) erste 
Paarungen bereits Ende Dezember. Der 
Hauptzeitpunkt ist aber Mitte Januar (DO¬ 
BOSZYNSKA & ZüROWSKl 1983). Im Untersu¬ 
chungsgebiet scheinen Beobachtungen von 
Jungtieren außerhalb des Baues Mitte Juni 
diesen Zeitpunkt für den Beginn der Paa¬ 
rungszeit ebenfalls zu bestätigen. 

Diese drei Faktoren - Zunahme der 
Markierhäufigkeit, Ende der Explorations¬ 
phase im Dezember des Aussetzungsjahres, 
Beginn der Paarungszeit - lassen es biolo¬ 
gisch sinnvoll und deshalb vertretbar er¬ 
scheinen, als Startmonat für ein „Biberjahr“ 
den Januar zu wählen und für Analysen im 
Zusammenhang mit dem Thema der Frei¬ 
landuntersuchung jeweils die Daten von Ja¬ 
nuar bis Dezember zusammenzufassen. 
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